Location: FTU Heart @ ad02dacd3e02 / Cardiac FTU v8.cellml
- Author:
- Soroush <ssaf006@aucklanduni.ac.nz>
- Date:
- 2021-10-01 19:36:57+13:00
- Desc:
- fixed the python script
- Permanent Source URI:
- https://models.cellml.org/workspace/70b/rawfile/ad02dacd3e022c6ffdcb7e80a6501170b07c2d75/Cardiac FTU v8.cellml
<?xml version='1.0' encoding='UTF-8'?>
<!--
Implements pressure calculation based on stress equilibrium but integrating from outside radius
Displays radial distributions of fibre angle and stresses and strains for radial expansion of an orthotropic thick-walled cylinder
Volume has units of L so radii are dm.
Pressure is J.L^-1;
2nd P-K stresses: T11 & T33 are J.L^-1; T22 is J.dm^-5; T23 is J.dm^-4 -->
<model name="Cardiac_FTU" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
<units name="dim">
<unit units="dimensionless"/>
</units>
<units name="dm">
<unit prefix="deci" units="metre"/>
</units>
<units name="dm2">
<unit exponent="2" prefix="deci" units="metre"/>
</units>
<units name="per_dm">
<unit exponent="-1" units="dm"/>
</units>
<units name="per_dm2">
<unit exponent="-2" units="dm"/>
</units>
<units name="radian_per_dm">
<unit units="radian"/>
<unit exponent="-1" units="dm"/>
</units>
<units name="J_per_L">
<unit units="joule"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="J_per_dm">
<unit units="joule"/>
<unit exponent="-1" units="dm"/>
</units>
<units name="J_per_dm2">
<unit units="joule"/>
<unit exponent="-2" units="dm"/>
</units>
<units name="J_per_dm4">
<unit units="joule"/>
<unit exponent="-4" units="dm"/>
</units>
<units name="J_per_dm5">
<unit units="joule"/>
<unit exponent="-5" units="dm"/>
</units>
<component name="main">
<!-- Integration variable-->
<variable initial_value="0" name="s" units="dimensionless"/>
<!-- Input Variables-->
<variable initial_value="0.0" name="T_0" units="J_per_L"/>
<variable initial_value="0.1" name="r_endo" units="dm"/>
<!-- Constants-->
<variable initial_value="0.1" name="R_endo" units="dm"/>
<variable initial_value="0.2" name="R_epi" units="dm"/>
<variable initial_value="0.0" name="phi_a" units="radian_per_dm"/>
<variable initial_value="1.0" name="lambda_a" units="dimensionless"/>
<variable name="lambda_f" units="dimensionless"/>
<variable initial_value="1" name="k_11" units="J_per_L"/>
<variable initial_value="1" name="k_22" units="J_per_dm"/>
<variable initial_value="1" name="k_33" units="J_per_L"/>
<variable initial_value="1" name="k_23" units="J_per_dm2"/>
<variable initial_value="1.2" name="c_11" units="dimensionless"/>
<variable initial_value="0.012" name="c_22" units="dm2"/>
<variable initial_value="1.2" name="c_33" units="dimensionless"/>
<variable initial_value="0.12" name="c_23" units="dm"/>
<!-- State Variables-->
<variable name="alpha_f" units="radian"/>
<variable name="r_epi" units="dm"/>
<variable name="r" units="dm"/>
<variable name="R" units="dm"/>
<variable name="E_rr_endo" units="dimensionless"/>
<variable name="E_rr_epi" units="dimensionless"/>
<variable name="E_11" units="dimensionless"/>
<variable name="E_22" units="dm2"/>
<variable name="E_33" units="dimensionless"/>
<variable name="E_23" units="dm"/>
<variable name="a11" units="dimensionless"/>
<variable name="a22" units="per_dm2"/>
<variable name="a33" units="dimensionless"/>
<variable name="a23" units="per_dm"/>
<variable name="Q_1" units="dimensionless"/>
<variable name="Q_2" units="J_per_L"/>
<variable name="Q_3" units="J_per_L"/>
<variable initial_value="0" name="f_23" units="dimensionless"/>
<variable initial_value="0" name="f_33" units="dimensionless"/>
<variable name="T11" units="J_per_L"/>
<variable name="T22" units="J_per_dm5"/>
<variable name="T33" units="J_per_L"/>
<variable name="T23" units="J_per_dm4"/>
<variable initial_value="0" name="p_0" units="J_per_L"/>
<variable name="p_epi" units="J_per_L"/>
<variable name="p" units="J_per_L"/>
<variable name="p_LV" units="J_per_L"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>r_epi</ci>
<apply>
<root/>
<apply>
<plus/>
<apply>
<power/>
<ci>r_endo</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<power/>
<ci>R_epi</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>R_endo</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<ci>lambda_a</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>R</ci>
<apply>
<plus/>
<ci>R_epi</ci>
<apply>
<times/>
<ci>s</ci>
<cn cellml:units="dim">0.001</cn>
<apply>
<minus/>
<ci>R_endo</ci>
<ci>R_epi</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>r</ci>
<apply>
<plus/>
<ci>r_epi</ci>
<apply>
<times/>
<ci>s</ci>
<cn cellml:units="dim">0.001</cn>
<apply>
<minus/>
<ci>r_endo</ci>
<ci>r_epi</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_rr_endo</ci>
<apply>
<times/>
<cn cellml:units="dim">0.5</cn>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<ci>R_endo</ci>
<apply>
<times/>
<ci>lambda_a</ci>
<ci>r_endo</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<cn cellml:units="dim">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_rr_epi</ci>
<apply>
<times/>
<cn cellml:units="dim">0.5</cn>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<ci>R_epi</ci>
<apply>
<times/>
<ci>lambda_a</ci>
<ci>r_epi</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<cn cellml:units="dim">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p_epi</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<ci>R_epi</ci>
<apply>
<times/>
<ci>lambda_a</ci>
<ci>r_epi</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<ci>k_11</ci>
<ci>E_rr_epi</ci>
</apply>
<apply>
<power/>
<apply>
<minus/>
<ci>c_11</ci>
<apply>
<abs/>
<ci>E_rr_epi</ci>
</apply>
</apply>
<cn cellml:units="dim">3</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_11</ci>
<apply>
<times/>
<cn cellml:units="dim">0.5</cn>
<apply>
<minus/>
<apply>
<power/>
<apply>
<divide/>
<ci>R</ci>
<apply>
<times/>
<ci>lambda_a</ci>
<ci>r</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<cn cellml:units="dim">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_22</ci>
<apply>
<times/>
<cn cellml:units="dim">0.5</cn>
<apply>
<minus/>
<apply>
<power/>
<ci>r</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>R</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_1</ci>
<apply>
<minus/>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<power/>
<apply>
<divide/>
<apply>
<times/>
<ci>r</ci>
<ci>phi_a</ci>
</apply>
<ci>lambda_a</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<apply>
<power/>
<ci>lambda_a</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dim">3</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>r</ci>
<ci>R</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dim">2</cn>
<ci>lambda_a</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>r</ci>
<ci>R</ci>
</apply>
<cn cellml:units="dim">4</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_2</ci>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<power/>
<ci>r</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>k_22</ci>
<ci>E_22</ci>
</apply>
<apply>
<power/>
<apply>
<minus/>
<ci>c_22</ci>
<apply>
<abs/>
<ci>E_22</ci>
</apply>
</apply>
<cn cellml:units="dim">3</cn>
</apply>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>k_11</ci>
<ci>E_11</ci>
</apply>
<apply>
<power/>
<apply>
<minus/>
<ci>c_11</ci>
<apply>
<abs/>
<ci>E_11</ci>
</apply>
</apply>
<cn cellml:units="dim">3</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_3</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<ci>k_11</ci>
<apply>
<power/>
<ci>lambda_a</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<plus/>
<ci>c_11</ci>
<apply>
<times/>
<cn cellml:units="dim">2</cn>
<ci>E_11</ci>
</apply>
</apply>
</apply>
<apply>
<power/>
<apply>
<minus/>
<ci>c_11</ci>
<apply>
<abs/>
<ci>E_11</ci>
</apply>
</apply>
<cn cellml:units="dim">4</cn>
</apply>
</apply>
<apply>
<minus/>
<ci>lambda_a</ci>
<apply>
<power/>
<apply>
<divide/>
<ci>R</ci>
<ci>r</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>s</ci>
</bvar>
<ci>p_0</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dim">0.001</cn>
<apply>
<minus/>
<ci>r_endo</ci>
<ci>r_epi</ci>
</apply>
</apply>
<ci>r</ci>
</apply>
<apply>
<power/>
<apply>
<divide/>
<ci>R</ci>
<apply>
<times/>
<ci>r</ci>
<ci>lambda_a</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>p</ci>
<ci>Q_1</ci>
</apply>
<ci>Q_2</ci>
</apply>
<ci>Q_3</ci>
</apply>
<apply>
<times/>
<ci>T_0</ci>
<apply>
<power/>
<apply>
<cos/>
<ci>alpha_f</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>p</ci>
<apply>
<plus/>
<ci>p_epi</ci>
<ci>p_0</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_f</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<minus/>
<pi/>
</apply>
<cn cellml:units="dim">3</cn>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>s</ci>
<cn cellml:units="dim">0.001</cn>
<cn cellml:units="dim">5</cn>
</apply>
<cn cellml:units="dim">6</cn>
</apply>
<pi/>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>lambda_f</ci>
<apply>
<root/>
<apply>
<plus/>
<apply>
<power/>
<apply>
<times/>
<ci>lambda_a</ci>
<apply>
<sin/>
<ci>alpha_f</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<times/>
<apply>
<power/>
<ci>r</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<cos/>
<ci>alpha_f</ci>
</apply>
<ci>R</ci>
</apply>
<apply>
<times/>
<ci>phi_a</ci>
<apply>
<sin/>
<ci>alpha_f</ci>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_33</ci>
<apply>
<times/>
<cn cellml:units="dim">0.5</cn>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<power/>
<ci>lambda_a</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<apply>
<times/>
<ci>r</ci>
<ci>phi_a</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<cn cellml:units="dim">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>E_23</ci>
<apply>
<times/>
<cn cellml:units="dim">0.5</cn>
<ci>phi_a</ci>
<apply>
<power/>
<ci>r</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>a11</ci>
<apply>
<power/>
<apply>
<divide/>
<apply>
<times/>
<ci>lambda_a</ci>
<ci>r</ci>
</apply>
<ci>R</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>a22</ci>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dim">1</cn>
<apply>
<power/>
<ci>r</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<power/>
<apply>
<divide/>
<ci>phi_a</ci>
<ci>lambda_a</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>a33</ci>
<apply>
<divide/>
<cn cellml:units="dim">1</cn>
<apply>
<power/>
<ci>lambda_a</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>a23</ci>
<apply>
<divide/>
<ci>phi_a</ci>
<apply>
<power/>
<ci>lambda_a</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>T11</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>k_11</ci>
<ci>E_11</ci>
</apply>
<apply>
<power/>
<apply>
<minus/>
<ci>c_11</ci>
<apply>
<abs/>
<ci>E_11</ci>
</apply>
</apply>
<cn cellml:units="dim">3</cn>
</apply>
</apply>
<apply>
<times/>
<ci>p</ci>
<ci>a11</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>T22</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>k_22</ci>
<ci>E_22</ci>
</apply>
<apply>
<power/>
<apply>
<minus/>
<ci>c_22</ci>
<apply>
<abs/>
<ci>E_22</ci>
</apply>
</apply>
<cn cellml:units="dim">3</cn>
</apply>
</apply>
<apply>
<times/>
<ci>p</ci>
<ci>a22</ci>
</apply>
<apply>
<times/>
<ci>T_0</ci>
<apply>
<power/>
<apply>
<cos/>
<ci>alpha_f</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>T33</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>k_33</ci>
<ci>E_33</ci>
</apply>
<apply>
<power/>
<apply>
<minus/>
<ci>c_33</ci>
<apply>
<abs/>
<ci>E_33</ci>
</apply>
</apply>
<cn cellml:units="dim">3</cn>
</apply>
</apply>
<apply>
<times/>
<ci>p</ci>
<ci>a33</ci>
</apply>
<apply>
<times/>
<ci>T_0</ci>
<apply>
<power/>
<apply>
<sin/>
<ci>alpha_f</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>T23</ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>k_23</ci>
<ci>E_23</ci>
</apply>
<apply>
<power/>
<apply>
<minus/>
<ci>c_23</ci>
<apply>
<abs/>
<ci>E_23</ci>
</apply>
</apply>
<cn cellml:units="dim">3</cn>
</apply>
</apply>
<apply>
<times/>
<ci>p</ci>
<ci>a23</ci>
</apply>
<apply>
<times/>
<ci>T_0</ci>
<apply>
<cos/>
<ci>alpha_f</ci>
</apply>
<apply>
<sin/>
<ci>alpha_f</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>s</ci>
</bvar>
<ci>f_33</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dim">0.002</cn>
<pi/>
<ci>r</ci>
<ci>T33</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>s</ci>
</bvar>
<ci>f_23</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dim">0.002</cn>
<pi/>
<ci>r</ci>
<ci>T23</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>p_LV</ci>
<ci>T11</ci>
</apply>
</math>
</component>
</model>